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Abstract. Multi-site and multi-year differential photometry of the triple star DG Leo reveals a complex frequency spectrum that 
can be modelled as the combination of at least three S Scuti type frequencies in the range 1 1.5-13 c/d (with semi-amplitudes of 
2-7 mmag) and a superimposed slow variability of larger amplitude. The period of the slow variation fits very well with half the 
orbital period of the inner spectroscopic binary indicating the presence of ellipsoidal variations caused by the tidally deformed 
components in a close configuration. These findings, together with the results of a recent spectroscopic analysis (showing that 
the system consists of a pair of mild Am stars and one A-type component of normal solar composition), infer that DG Leo is 
an extremely interesting asteroseismic target. Identification of which component(s) of this multiple system is (or are) pulsating 
and determination of the excited pulsation modes will both contribute to a much better understanding of the non-trivial link 
between multiplicity, chemical composition, rotation, and pulsation in the lower part of the classical Cepheid instability strip. 
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Technique: photometry 



1. Introduction 

1.1. Objectives 

A 'bonus' of studying stellar pulsations in binary and multi- 
ple systems is that one can exploit the dynamical information 
to obtain independent constraints on the pulsational models via 
the determination of the physical properties of the components. 
Multi-year photometry with a good phase coverage of the beat 
periods of the oscillations, on the one hand, allows accurate 
frequency analyses of the combined light. A careful spectro- 
scopic analysis will, on the other hand, provide information on 
the mass and the luminosity ratios, the effective temperature, 
the chemical composition, and the surface gravity of each com- 
ponent given that the spectrum is composite. If dynamically de- 
termined masses of the system can be obtained with sufficiently 
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* Based on photoelectric observations obtained at San Pedro Martfr, 
Sierra Nevada observatories, and CCD observations at Beersel Hills 
Observatory, as well as on photoelectric observations with the Belgian 
mercator telescope 



good accuracy, these can be directly confronted with the pulsa- 
tion masses. By comparing the frequency and modal content of 
the components that originate from the same protostellar envi- 
ronment, one may realistically hope to improve understanding 
of pulsation physics, since a difference in pulsational behaviour 
between each can only be attributed to a limited number of (dif- 
fering) stellar parameters. Binarity/multiplicity may thus lead 
to an explanation of the so-far unknown amplitude and mode 
selection mechanism operating in 5 Scuti stars, and to a refined 
characterization of the pulsations of variable stars in general. 

In a review of pulsating 6 Scuti stars in known double 
or multiple stars, Lampens & Boffin (2001) discussed several 
potentially interesting objects, one of which was DG Leonis 
which appears to be a triple system whose components are 
nearly identical. 

1.2. The target 

DG Leo (HR 3889, HIP 48218, Kui 44 AB) is an hierarchical 
triple system which consists of a close double-lined spectro- 
scopic binary (components Aa and Ab) and one distant com- 
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panion (component B) forming the wider visual binary. It has V 
= 6.085 mag and (B-V) = 0.25 mag (Mermilliod & Mermilliod 
1998 1. The orbital period of the Aa,b system is 4.147 days, 
while the orbital period of the visual pair AB is roughly 200 
years (Fekel & Bopp 1977). The visual component B, ob- 
served by hipparcos with an angular separation of 0.17"and a 
magnitude difference of 0.69 mag with respect to Aa,b (ESA 
1997), was monitored by means of micrometric observations 
and speckle interferometry from 1935 till 1997 (Hartkopf et 
al. |2000l. The existing astrometric data indicate a highly in- 
clined and possibly eccentric orbit. Both Barnes et al ([1977 1 
and Rosvick & Scarfe (1991 1 presented evidence for possible 
shallow eclipses with an inclination of w 75° for the Aa,b sys- 
tem. 

Both close components have spectral type A8 IV (Hoffleit 
& Jaschek ll982> . while the composite spectrum was classified 
as type F0 Illn (Cowley 1976 1 and type A7 III with enhanced Sr 
(Cowley & Bidelman 1979 1. All three components of the sys- 
tem are located in the 6 Scuti instability strip and are therefore 
potential candidates for pulsations. Danziger & Dickens ( 1967 1 
first reported the short-period variability in the system but were 
inconclusive about its nature. The claims for short-period os- 
cillations seem to concern mainly one component: component 
B was classified both as an ultra-short period Cepheid (Eggen 
119791 with a period < 0.275 days) and as a 6 Scuti star (Elliot 
119741 with a mean periodicity of 0.0818 days). Possible am- 
plitude and phase changes were also reported (Antonello & 
Mantegazza 1982 Rosvick & Scarfe 1991 1. Since these pre- 
vious studies are based on very limited photometric data sets, 
little was known about the pulsational characteristics. A first 
report on the multiperiodicity of the variations was recently 
presented after performing a preliminary frequency analysis 
of newly acquired photometric data (Lampens et al. 2004). In 
the following sections, we will present the complete and up- 
dated photometric analysis of this interesting asteroseismic tar- 
get, where metallicity as well as binarity effects may influence 
pulsational behaviour. 

2. Observations and data reduction 

2.1 . Photoelectric photometry 

During 2002 and 2003, we performed photoelectric photom- 
etry campaigns of DG Leo at two prime sites located at dif- 
ferent longitudes and equipped with identical instrumentation 
and standard Stromgren filters in a near simultaneous mode. 
New uvby data were collected with the 1.5 m telescopes at San 
Pedro Martfr (OAN), Mexico, between Feb. 4 - 11, 2002 and 
at Sierra Nevada (OSN), Spain, between Feb. 11 - 17, 2002. 
More uvby data were collected at OAN between January 8 - 14, 
2003. We used CI = HD 84497 = HIP 47946 = GSC 1416-210 
(G8 III, V=7.42, b-y=0.581) and C2 = HD 84739 = HIP 4805 1 
= Cou 284 (F0, V=7.39, b-y=0.290, also a close visual double 
star) as comparison stars. The check star was C3 = HD 86516 
= HIP 48977 = GSC 1418-1085 (Am, V=6.74, b-y=0.109). 
The magnitude differences in the sense (CI - C3) were shown 
by Rosvick & Scarfe (1991) to be constant at their level of 
precision. The hipparcos catalogue (ESA I1997t mentions that 



CI showed no detectable variability throughout the mission, 
while the comparison stars C2 and C3 carry the flag 'Duplicity- 
induced variability'. The standard deviations of the uvby data 
of the OSN 2002 and OAN 2003 campaigns show that CI 
and C2 remained constant at the level of a few mmag, ex- 
cept in filter u for which the noise level is larger (cf. values 
between brackets in Table 2 introduced later in this Sect.). On 
the other hand, C3 might be variable on a long-term scale, even 
though it was used for a short time as a comparison star by 
Jerzykiewicz ( 1993 1. We found a difference of more than 2 <x 
on index ml and more than 6 <x on index cl when compar- 
ing our mean indices with the standard indices published by 
Maitzen et al. ( 1998, based on 6 measurements, see below). 

The observations at OAN were performed according to the 
sequence 'variable sky Cl variable C2 ( variable C3 )'. Five 
successive 10-second integrations per filter were used for the 
stars, while only one 10-second integration per filter was used 
for the sky. The observations at OSN were performed accord- 
ing to the sequence 'variable Cl C2 sky' using 30-second in- 
tegrations in each filter for all objects. During the OAN 2003 
campaign the check star (C3) was also included in the observa- 
tional sequence. 

All OAN data were reduced according to the same reduc- 
tion method as applied to the OSN data. First, the dead-time 
correction was applied to all counts. Then, subtraction of the 
(average) sky measurement from the average value of the tar- 
get star was carried out for each cycle. The atmospheric extinc- 
tion was calculated from the comparison star data following 
the method of Gr0nbech et al. (1976) and the corresponding 
corrections were subsequently applied to the magnitude differ- 
ences. The OAN instrumental errors computed for each season 
using the standard stars are shown in Table 1. As can be seen 
from this table, the OAN instrumental errors were larger for 
the 2002 than for the 2003 season. The reason for the smaller 
errors in 2003 was the implementation of new phototubes for 
the uvby photometric system in 2002 followed by a debugging 
phase which lead to more precise determination of the dead- 
time corrections for the 2003 OAN season. 

The photometric data discussed here represent magnitude 
differences in the instrumental system between the variable and 
the comparison star (Cl or C2) interpolated to the observa- 
tion time of the variable star. As the mean differences were 
not exactly the same for both observatories, we applied small 
corrective shifts for each of the three campaigns (see Table 2). 
Note that the OAN 2003 campaign was held at an epoch closely 
matching that of the the high-resolution spectroscopic run con- 
ducted at the Observatoire de Haute-Provence (Fremat et al. 
2004b I. In this way we obtained a homogeneous data set con- 
sisting of 1590 epochs and multi-colour differential magnitudes 
of DG Leo. The new material represents about 110 hrs of high- 
quality photoelectric photometry (Data Set 1). 

2.2. CCD photometry 

In order to increase the amount of data and extend the time 
basis of the previous data set, we also performed comple- 
mentary observations with a small telescope at Beersel Hills 
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Dates 


Observatory 


Observer(s) 


Nr 


Comp 1 


Comp 2 


Check 


His 


Feb. 4- 11,2002 


OAN 


JP 


429 


CI 


C2 


no 


40.6 


Feb. 11 - 17, 2002 


OSN 


RG 


748 


CI 


C2 


no 


42.0 


Feb. 22 -Mar. 18, 2002 


LPA 


PDC & KU 


90 


CI 


no 


no 




Jan. 8- 14, 2003 


OAN 


JP 


269 


CI 


C2 


C3 


20.0 


Jan. 9 - Mar. 14, 2003 


BHO 


PL & PVC 


1071 


C4 


no 


C5 


58.7 


Mar. 17- 19, 2004 


OSN 


RG 


144 


CI 


C2 


no 


7.0 



*: non-continuous observations obtained at La Palma observatory (LPA), cf. Sect. 4 



Observatory (BHO), Belgium, during the year 2003. Since 
the target is very bright, we used a 10cm telescope with an 
SBIG STIOXMe camera to additionally collect differential 
CCD data, and we adopted C4 = HD 85017 = HIP 48197 = 
GSC 1417-315 as the comparison star and C5 = HD 85215 
= HIP 48303 = GSC 1417-312 as the check star. We also 
defocused slightly in order to improve the photometric preci- 
sion. These frames were reduced with the aperture photometry 
procedure of the Mira AP software package 1 . 1071 differential 
(stacked) V measurements were thus obtained between Jan. 9 
and March 14, 2003 representing a total of 59 hours of differ- 
ential CCD photometry in the instrumental system. Though the 
equipment and the filters are not exactly the same and the com- 
parison stars were obviously different, we combined all dif- 
ferential V and Stromgren y measurements into one large set 
(Data Set 2 containing 2661 differential magnitudes in a 'sin- 
gle' filter 2 ), which will serve the purposes of a) providing a dif- 
ferent time distribution with a distinct alias pattern and b) ver- 
ifying the frequency solution that will be presented in the fol- 
lowing section. To perform such data combination, we needed 
to apply small corrective shifts of order 0.01 mag on the CCD 
differential data computed from the shifts between the nightly 
means of the two comparison stars in the field. The averaged 
V and Stromgren y magnitude differences were subsequently 
subtracted for both data types. No other transformation was 
applied. Table 3 shows the log of the observational runs per- 
formed between January 2002 and March 2004. 

3. Analyses 

3.1. Frequency analysis of Data Set 1 

We performed the frequency analysis with two different codes 
that are based on computation of the Fourier Tranform and/or 
on the technique of least-squares fitting. 

At first, we used Period98 (Sperl 1998), a programme 
which determines the frequencies one by one by computing 
the classical Fourier Transform followed by adjustment of the 
parameters of the periodic functions (zero point, amplitudes, 
and phases) using least-squares fitting. Next, the computed fre- 
quency solution is refined through simultaneous adjustment of 

1 The software Mira AP is produced by Axiom Research Inc., 
http://www.axres.com/ 

2 The data will be available in electronic form via anony- 
mous ftpto cdsarc.u-strasbg.fr (130.79.128.5) or via 
http : //cdsweb . u-strasbg . fr/cgi-bin/qcat? J/A+A/xxx/xxx 



Table 1. OAN errors in the instrumental system and colour in- 
dices of the comparison stars in the standard system (season 
2003) 



Season 


Filter u 


Filter v 


Filter b 


Filter y 


OAN 2002 


0.031 


0.012 


0.008 


0.008 


OAN 2003 


0.024 


0.005 


0.006 


0.006 


Star 


V(oy) 


(b-y) (ov,) 


ml (<r m i) 


cl (cr c i) 


CI 


7.394 (20) 


0.583 (11) 


0.326 (28) 


0.346 (45) 


C2 


7.358 (20) 


0.291 (11) 


0.144 (28) 


0.473 (45) 


C3 


6.729 (20) 


0.097(11) 


0.282 (28) 


0.563 (45) 



all the parameters, including the frequencies in the course of a 
(local) minimization process. Since we are dealing with mul- 
tiperiodic signals, we used the empirical criterion that a fre- 
quency is significant as long as the signal-to-noise ratio of the 
corresponding semi-amplitude is greater than 4 (Kuschnig et 
al. lTWTl . 

Secondly, a multi-parameter fitting code that scans a wide 
range in frequency and searches to minimize the residuals by 
simultaneously adjusting multiple periodic functions was used 
(global minimization). In this way accurate, alternative multi- 
frequency solutions can be found that might be overlooked by 
the previous code (Schoenaers & Cuypers 2004 1. This is espe- 
cially useful when an alternative solution due to the ambigu- 
ity of some aliasing pattern is also probable. In this case both 
computations gave identical results in all filters except for filter 
u which was of lower quality. Table 4 lists the results for Data 
Set 1 for each filter of the Stromgren system. Listed are the 
frequency (in cd~'), the semi-amplitude (in mmag), the phase 
(in 27rrad), and the total fraction of the variance that is removed 
from the signal (R). In addition we list the signal-to-noise ratio, 
with respect to the mean noise level in the residual periodogram 
computed within a box of width 5c d _1 centered on the relevant 
frequency, as well as the standard deviation of the residual sig- 
nal, cr res , which would correspond if Period98 were used; and a 
succession of prewhitenings with the previous frequency would 
be performed at each step. 

The dominant frequency Fl, the one of largest amplitude 
found at about 0.48 cd~', corresponds to a slow periodic 
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change of the mean intensity linked to the orbital variations 
of the short-period pair Aa, b. It is very close to one half of 
Porb - 4.146751 days, which is the orbital period derived from 
the published spectroscopic analyses (Fekel & Bopp 119771 
Fremat et al. 2004a). Note that the orbital frequency itself (F0) 
is not significant, as its signal-to-noise ratio remains far below 
4 in all the passbands. However, we considered its contribu- 
tion, since omitting it would affect the other frequency values 
very slightly, particularly in the u passband. Four more frequen- 
cies between 1 1.9 and 12.8 c d~'(F2-F5) were detected in each 
filter. In fact, two of these are located about 0.08 cd~' apart, 
which caused some caution in the preliminary study where T 
was only about 12 days (Lampens et al. 2004 1. Their presence 
is undoubtedly confirmed as the new time base of 767 days cor- 
responds to a formal resolution of 0.0013 cd -1 . Also note that 
F3 and F4 could be subject to the n.year -1 ambiguity (with n 
equal to 1 or 3). Obviously, there is excellent agreement be- 
tween the solutions in all four passbands. We point out the high 
R-value in filters b and v, but also the poorer match of the same 
solution in the u passband (R ~ 33%). From a comparison of 
the R-values for nearby alias frequency sets, we verified that no 
other combination was as significant as the one listed for the y 
passband. As an illustration of this analysis we present Fig. [2] 
which is entirely based on the computations performed in par- 
allel with Period98 and derived in successive steps for the data 
in filter b. 

The expected error in frequency in the case of uncorrelated 
observations can be computed using the equation 



a F 



V6 



n 'A.N l ' 2 .T 



(Cuypers 1987 1, where A/cr„ indicates the signal-to-noise ra- 
tio, N, the number of observations and T, the total time base. 
From Table 4, we adopt 5 mmag as the typical noise level re- 
maining after prewhitening for six frequencies in the filters v, 
b, and y and derive an error in frequency of l.-5.10~ 5 cd _1 on 
frequencies F1-F5, respectively, given that T = 767 days and 
N = 1590 (Data Set 1). This also corresponds to the error one 
would obtain adopting the formula derived by Schwarzenberg- 
Czerny (I199U . provided that some of the simplifications are 
not considered. As noted by Schwarzenberg-Czerny ( 1991 1 and 
Montgomery & O'Donoghue (1999 1, correlations in the noise 
may increase this error with the square root of the noise correla- 
tion length V7J. To estimate D, we analysed the residuals after 
full prewhitening with six frequencies and found a mean sign 
change at almost every 4 data points. In the following we will 
adopt the conservative error estimates of 2.-10. 10 -5 c d'on fre- 
quencies F1-F5. In view of the small differences obtained be- 
tween the frequency analyses in the filters y, b, and v, these 
errors further suggest that a unique solution is valid for all the 
data. Apart from these values, which can be considered as up- 
per limits of the errors on the frequencies, we recall the fact 
that the n.year -1 aliasing effect still affects the current analysis 
and that this is a practical limitation to the true accuracy of the 
frequencies F3 and F4 in Table 4. We decided to stop the fre- 
quency search at this level since we noted a significant change 
of the slope at N = 4 when plotting the residual standard de- 



Fig. 1. Residual standard deviation as a function of the number 
of removed sinusoids (Set 1 - filter b) 

Table 4. Multi -parameter fits and solutions (Data Set 1) 



Set 


Freq. 


Amp. 


Phase 




S/N 


& res 


R 


1 


ccT 1 


mmag 


27rrad 






mmag 


% 




Filter y 


N 


= 1590 


0~init = 


8.4 mmag 






Fl 


0.482282 


7.7 


0.41 




22 


0.7 




(F0 


0.241141 


0.8 


0.36 




2.2 


- \ 

i 




F2 


11.99429 


5.2 


0.98 




13 


5.8 




F3 


11.91064 


4.5 


0.63 




1 1 


5.0 




F4 


12.11436 


3.2 


0.00 




8 


4.5 




F5 


12.77599 


1.6 


0.32 




4.0 


4.3 


74 




Filter b 


N 


=1590 


0~inir = 


9.5 mmag 






Fl 


0.482276 


8.4 


0.41 




15 


7.7 




(F0 


0.241138 


0.4 


0.31 




0.8 


-) 




F2 


11.99426 


6.2 


0.97 




13 


6.5 




F3 


11.91062* 


5.3 


0.62 




11 


5.4 




F4 


12.11439" 


3.8 


0.02 




8 


4.7 




F5 


12.77592 


2.3 


0.34 




4.9 


4.4 


79 




Filter v 


N 


= 1590 


& intt~ 


10.8 mmag 






Fl 


0.482274 


8.9 


0.41 




14 


9.0 




(F0 


0.241137 


0.4 


0.38 




0.6 


-) 




F2 


11.99427 


7.3 


0.97 




13 


7.6 




F3 


11.91062* 


6.4 


0.62 




12 


6.2 




F4 


12.11439" 


4.7 


0.02 




9 


5.3 




F5 


12.77591 


2.6 


0.34 




4.8 


4.9 


79 




Filter u 


N 


= 1590 


C imt = 


15.5 mmag 






Fl 


0.482268 


8.1 


0.41 




6 


14.3 




(F0 


0.241134 


2.8 


0.57 




2.1 


-) 




F2 


11.99425 


6.5 


0.95 




7 


13.7 




F3 


11.91055 


6.3 


0.59 




7 


13.1 




F4 


12.11432 


3.7 


0.94 




4 


12.8 




F5 


12.77586 


2.9 


0.35 




3.2 


12.7 


33 



* alternatively, the 3. (year) 1 alias frequency 11.9197 cd 1 
** alternatively, the l.(year)~' alias frequency 12.1172 cd" 1 



viation of the b data as a function of the number of sinusoids 
removed from the original data (Fig.0. 

3.2. Frequency analysis of Data Set 2 

Data Set 2 spans the same time base as the former data set 
but consists of 2661 combined measurements in the filters y 
and V. While the time base is identical, the increased number 




0.008 E- 
0.006 E- 
0.004 - 



i i i i l i i i i l i i i i l i i i i l i i i i . 



I 



£ 0.006 E- 

^ 0.004 E- 
1) 

0.002 

> 

5 

5 

flj r-r 

^ 0.008 F 
g 0.006 

6 0.004 
< 0.002 





10 15 20 25 

0.008 I 1 
0.006 E- 
0.004 E- 
0.002 



0X02 taittiM^iiirtJtfifliA^^l^A^ 



F5 



10 15 20 25 



i i i i I i i i i I i i i i I i i i i I i i i i . 





10 15 20 25 



5 10 15 20 25 




1 

0.8 



0.2 




J 1 II II 1 II 


1 1 1 1 1 1 1 1 1 1 




Spectral window ~ 


r ,A 





-10 







10 



I I I I I I I I I I I I I I I I I I I I I I I I I [ 

Spectral window ~ 



10 15 20 25 




-0.1 -0.05 0.05 0.1 



Frequency (c/d) 

Fig. 2. Successive frequency searches and spectral window of the Fourier analysis (Set 1 - filter b) 



of measurements should improve the accuracy of the analysis. 
However, because the site for the CCD data is located near an 
industrialized area and a smaller telescope was used, the accu- 
racy of a single measurement is lower. We therefore assigned a 
relative weight of about 5 _1 for each CCD measurement with 
respect to a photoelectric one. The inverse ratio of the variances 
of the differential magnitudes between the comparison stars C4 
and C5 (BHO) and CI and C2 (OSN, OAN) was used to com- 
pute this factor. The main asset of these data is that the spectral 
window is different, largely due to the fact that the CCD obser- 
vations spanned several months in the 2003 season (cf. Fig.0- 

We again performed the same analysis twice: first with 
Period98, then with the method of multi-parameter fitting. We 
list the best set of parameters when adopting a six-frequency 
model with the latter method in Table 5. From a comparison of 
the R-values for close sets of aliased frequencies, we also veri- 
fied that no other combination was as significant as this one. In 
particular, the frequency at 12.1145 cd~'(F4) is slightly more 
convincing than its year -1 alias at 12. 1 172 c d _1 . Note that the 
3.year~' ambiguity no longer affects the solution. In theory, 
the formal errors in frequency would be improved by y/Ni /N2 
(with Ni and N2 the number of data in the first and last set re- 
spectively). However, considering the use of a relative weight 
of 5 _I for the V data, we will adopt the same conservative er- 



Table 5. Multi-parameter fit and solution (weighted Data Set 2) 



Set 


Freq. Amp. 


Phase 


S/N 


Vres 


R 


2 


cd _1 


mmag 






mmag 


% 




Filter y+V 


N 


=2661 


cr,„„= 9. 1 






Fl 


0.482278 


7.6 


0.42 


13 


7.5 




(F0 


0.241139 


1.1 


0.23 


2.7 


-) 




F2 


11.99425 


5.0 


0.98 


13 


6.7 




F3 


11.91064 


4.0 


0.63 


10 


6.1 




F4 


12.11446* 


2.9 


0.00 


8 


5.8 




F5 


12.77588 


1.6 


0.31 


4.2 


5.6 


62 



* alternatively, the l.(year) 1 alias frequency 12.1172 cd 1 



rors as before (i.e. errors of 2.-10. 10 -5 c/d respectively on fre- 
quencies F1-F5). From here on, we will adopt the frequency 
solution of the v data (Data Set 1) as the final solution since it 
is based on homogeneous data of very high quality and since 
it is furthermore fully consistent with the solution derived for 
Data Set 2 within the quoted error bars. 




0.1 -0.1 
Frequency (c/d) 



1 I . 

Spectral window 2 




Fig. 3. Spectral windows showing the difference in aliasing between Set 1 and Set 2 



3.3. Interpretation of the light variability 

As mentioned, the frequency of highest amplitude has a value 
close to 2.f m -h and corresponds to a slow periodic variability 
of the daily mean intensity due to the geometrical variations of 
the inner binary Aa, b. With an orbital period as short as 4 days, 
both components of Aa, b must have their shapes distorted by 
tidal interaction. We adopted the spectroscopic orbital period 
(?orb - 4.146751 days) and used the following model of the 
light curve of an ellipsoidal binary: 

LaoA'P) = a +/5 * cos((p) + y * cos(2<p) (1) 

with <p the orbital phase and a (the mean light level), /3, y the 
coefficients of the fit (Beech [T989> . The term corresponding 
to 2 * fori, (Fl) can be associated to the ellipticity effects: the 
non-spherical cross-sections vary as the system rotates around 
the centre of mass and maximum luminosity occurs twice per 
revolution. The term corresponding to / or j (F0) is related to 
the irradiation effect which occurs once per revolution: the in- 
ner part of the cooler component is irradiated by the hotter one 
and re-emits the excess radiation isotropic ally. There is no net 
observable effect in this case, meaning that both components 
must have similar effective temperatures (Wilson 1994). The 
question that remains is whether eclipses also occur or not. 
When fitting the OSN-OAN 2002 data with a single sinusoidal 
term as a function of half the orbital period, we note the al- 
most constant differences in the nightly means of the OAN run 
which can be assigned either to shallow eclipses and/or to el- 
lipsoidal variations (Fig. 4), whereas the slow change of the 
nightly mean levels in the OSN run is evidence for ellipticity 
effects (Fig. 5). This fit also shows that no (shallow) eclipses are 
needed to explain the observed light variations, in contrast to a 
previous estimation (Rosvick & Scarfe [T991> . Thus, the dom- 
inant photometric variation can be explained as caused solely 
by the tidally deformed shapes of both components. 

The other frequencies detected in the light variations may 
be associated to pulsations of the type 6 Scuti in one (possibly 
more) component(s) of DG Leo. Adopting a unique frequency 
solution (the final one derived in Sect. 3.2) for all passbands, 
we recomputed the semi-amplitudes and the phases with their 
standard errors through least squares fitting (Table 6). Note 
the tiny colour effect in frequencies Fl and F0 as the ampli- 
tude in filter u for F0 = f or b becomes larger. The amplitude 
ratios and phase shifts between the different colours displayed 
in this table contain information that will be helpful for identi- 
fying the modes of the excited pulsations (Garrido 2000). Note, 




Q 0.03 - 



JD - 2452317.042 

Fig. 4. Individual OAN light curves showing the ellipsoidal 
variation of DG Leo (filter y) 




Q 0.03 - 



JD - 2452317.042 

Fig. 5. Individual OSN light curves showing the ellipsoidal 
variation of DG Leo (filter y). Note the different orbital phase 
coverage. 

however, that the true semi-amplitudes will be larger than ob- 
served since three almost equally luminous stars were observed 
jointly. Given that A lr is the true semi-amplitude, A t, s the ob- 
served one, and f the additional light source with respect to the 
maximum light of the variable component B, we derive: 

1 + /.10 0M " 

A lr - A obs = 2.5 * log jj— (2) 

Application of this equation with f = 2 gives the corrected semi- 
amplitudes 0.021 and 0.006 mag corresponding to the extreme 
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Table 6. Comparison of amplitudes and phases in various pass- 
bands 



Set Filter u Filter v Filter b Filter y 



1 


A, 


©1 


A 2 


i>2 


A 3 


0>3 


A 4 


$4 




mmag 


2nraA 


mmag 


2;rrad 


mmag 


2;nad 


mmag 


2;rrad 


Fl 


8.1 


0.34 


8.9 


0.34 


8.4 


0.34 


7.6 


0.34 




±0.5 


0.06 


0.2 


0.02 


0.2 


0.02 


0.2 


0.02 


(FO 


2.8 


0.43 


0.4 


0.62 


0.4 


0.69 


0.8 


0.64) 




±0.5 


0.16 


0.2 


0.43 


0.2 


0.40 


0.2 


0.20 


F2 


6.3 


0.45 


7.3 


0.43 


6.2 


0.44 


5.3 


0.43 




±0.5 


0.07 


0.2 


0.02 


0.2 


0.03 


0.2 


0.03 


F3 


6.2 


0.20 


6.4 


0.17 


5.3 


0.17 


4.4 


0.16 




±0.5 


0.07 


0.2 


0.03 


0.2 


0.03 


0.2 


0.03 


F4 


3.7 


0.33 


4.7 


0.25 


3.8 


0.26 


3.1 


0.28 




±0.5 


0.12 


0.2 


0.04 


0.2 


0.04 


0.2 


0.05 


F5 


2.9 


0.72 


2.6 


0.74 


2.3 


0.75 


1.6 


0.77 




±0.5 


0.16 


0.2 


0.07 


0.2 


0.07 


0.2 


0.10 



values of 0.007 and 0.002 mag for the observed amplitudes of 
the pulsation, i.e. a factor of 3 larger. On the other hand, the 
semi-amplitudes related to the ellipsoidal variations should be 
corrected by an additional 50% (applying the same equation 
with f = 0.5 since the additional light source is component B 
with respect to components Aa,b). 

4. Determination of global atmospheric properties 

We next computed the OAN mean magnitudes and colours of 
the triple star in the Stromgren standard system (for the 2002 
season) in order to derive the global atmospheric parameters 
of the system through application of the calibration procedure 
(Lester, Gray & Kurucz 1986 1. The absolute magnitude, effec- 
tive temperature, surface gravity, and metallicity of DG Leo 
Aab,B and their standard errors are presented in Table 7. 

We also started to collect absolute data in all seven fil- 
ters of the Geneva photometric system with the mercator tele- 
scope 3 . We determined the mean colours and standard devia- 
tions of the 90 measurements acquired between February and 
March 2002 (cf. Table 1). Intrinsic colours of B2 to M0 stars 
in the Geneva system have been estimated by Meylan et al. 
(1980) and Hauck (1993). The global colours of DG Leo seem 
to match an unreddened MK spectral type of A7 III, except 
for filter U where the system is less luminous (cf. Table 2 in 
Hauck 1993 ). We next made use of the calib code developed 
by Kiinzli et al. (119971 1 which allows the atmospheric parame- 
ters from the Geneva colour indices to be derived for B to mid- 
G stars on or just above the main sequence. With this code it is 
also possible to treat metallic-line stars. Using a zero colour ex- 
cess, we obtained the atmospheric properties listed in Table 7. 
From the reddening-free parameter d = 1.283 (measuring the 
Balmer discontinuity), we next derived Ad = 0.107 (measuring 
how much the star is evolved away from the main sequence) 
and My = 1.51 ± 0.30 mag following Hauck ( 1973 1. However, 

3 The mercator telescope is operated by the Institute of Astronomy, 
K.U.Leuven, Belgium, at La Palma (Spain). 



North et al. ( 1997 1 re-examined the existing calibrations of the 
Stromgren and Geneva systems in terms of My for Am stars. 
Their conclusion was that the previous calibrations were not 
reliable and that the Geneva photometry underestimates the lu- 
minosity of the hottest Am stars. Applying the new relation on 
the Stromgren indices from Table 7 and using V sin i m 30 km/s 
as an appropriate value (Fremat et al. 2004b I, we derived the 
somewhat lower absolute magnitude My = 1.08 + 0.31 mag. 

We might compare these calibrated mean values with those 
that can be directly derived on the basis of the Hipparcos paral- 
lax (ESA 1997 1. In Sect. l3.3l we showed that both components 
of the close binary, Aa and Ab, have similar colours since no 
net reflection effect is observed in the light curve. We therefore 
assume that the difference in magnitude between these compo- 
nents is very small (neither one is evolved). On the other hand, 
we know that the (Hp) magnitude difference between compo- 
nents Aab and B is 0.69 mag (ESA [TW7l . Adopting V = 6.085 
mag, we then obtain 'SI Aab - 6.546 mag and = 7.236 mag. 
Similarly we also have Va = "V ' m> = 7.299 mag. Converting to 
absolute visual magnitudes gives My All = = 1.31 mag and 
My B = 1 .25 mag. The relative parallax error of 15% contributes 
the most to the error in such a way that none of the above cali- 
brations can be rejected (cr(My) w 2.17^ = 0.33 mag). 

The listed errors were increased by 50% with respect to the 
nominal accuracies to take the fact into account that three stars 
were observed jointly. Since the three components were not 
very different, these derived properties will correspond to ap- 
proximate values of the physical parameters of the components, 
whereas the high-resolution spectroscopic analysis allowed us 
to perform a detailed study of the differences between each of 
the components (cf. Table 6 in Fremat et al. 2004b). A compar- 
ison between these component properties and the global val- 
ues based on photometric calibrations shows that the agreement 
in effective temperature is reasonably good (mean of 7500 K), 
while the agreement in surface gravity is concordant with the 
value derived from the Geneva photometry (mean of 3.8). Both 
[M/H] determinations suggest enhanced metallicity of the com- 
bined light, an observational fact already mentioned by Fekel 
& Bopp (. 19771 and Cowley & Bildelman i 197% , but only re- 
cently quantified (Fremat et al. 2004b). 

Acquisition of absolute multi-colour photometry in the 
Geneva photometric system was pursued during 2004 with the 
aim of obtaining well-sampled colour curves for each of the 
pulsation frequencies. The reduction of these observations is 
presently taking place. 

5. Conclusions and future work 

The multi-colour time series of DG Leo recently collected on 
two continents at the same epoch were merged together and 
frequency-analysed. Both time series contain high-quality data 
in the Stromgren instrumental system. Their combination fur- 
thermore presents the advantage of largely suppressed day -1 
aliasing effects which made the frequency-analyses straightfor- 
ward and lead (a) to an improved interpretation of (the light 
curve of) the close binary system and (b) to the detection of 
multiperiodic oscillations in the combined flux. Monitoring 
during a few successive years (2002-2004) procured reasonably 
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good accuracy and also permitted verification of the stability of 
the new frequency solution. 

Up to four 6 Scuti frequencies in the range 11.5-13 cd _1 
(with semi-amplitudes of 2-7 mmag in the v passband), as well 
as a slow variation (with a semi-amplitude of 9 mmag in the v 
passband) were identified with high confidence. The frequency 
solution as derived from the Stromgren v data appeared to fit 
the data very well in all other passbands, as well as the more 
extended data set which is, however, less homogeneous. This 
unique frequency solution was shown to remain stable over a 
period of a few years. We further concluded that the period of 
the slow variation fits very well with half the spectroscopic or- 
bital period of the close binary and that its presence can be ex- 
plained by continuous ellipsoidal variations supplemented by 
a tiny colour effect. Despite these efforts, the full frequency 
content of the pulsations in DG Leo has not yet been identi- 
fied as still more peaks are present in the same range of the 
periodogram of the residuals. A far more extended time series 
would be needed to determine them with good reliability. 

As the triple system cannot be resolved photometrically, 
this analysis is valid for the combined light of all three com- 
ponents. Previous studies indicate that the components have 
similar characteristics. Therefore, the true semi-amplitudes of 
the pulsations are diluted by a factor of about 3. The semi- 
amplitudes and the phases were determined in each of the four 
Stromgren colour indices and will - after finalising the reduc- 
tion of the Mercator data - also be known in the Geneva photo- 
metric system. These will be important diagnostics when iden- 
tifying the degree t of the excited pulsation modes. Of course, 
such analysis depends on correct identification of the compo- 
nents) in which the pulsations originate. Only a detailed spec- 
troscopic analysis can provide reliable and accurate informa- 
tion at the component level. Fremat et al. ( 2004b I thus revealed 
that the close binary consists of two mild Am stars and that 
line-profile variations on the photometric time-scales of the 
pulsations were detected in the time-series spectra of the distant 
companion. The in-depth analysis of DG Leo - a triple system 
with three almost identical components - should put strict em- 
pirical constraints on the development of this type of pulsations 
in stars of the same mass and age, but of different chemical 
composition in the outer layers. 

Our intention is to postpone the modal interpretation until 
the frequency analysis of the line profile variations detected in 
the high-resolution spectra can be presented as well. Reliable 
conclusions on this issue should indeed come from evidence 
based on both photometric and spectroscopic diagnostics (both 
types of data were acquired quasi-simultaneously for a better 
understanding of this triple system). Only then will we be able 
to consider the exact origin of the pulsations and to establish the 
link between multiplicity, chemical composition, and pulsation 
in this extremely interesting S Scuti star. A modelling of the 
pulsations should finally be attempted, in the same way as the 
preliminary one by Grigahcene et al. (2005). 

Acknowledgements. P. Lampens, Y. Fremat, and T. Arentoft grate- 
fully acknowledge financial support from the Belgian Science Policy 
(through project MO/33/007). R. Garrido aknowledges financial sup- 
port from the project ESP2001-4528-EP. L. Parrao and J.H. Pena 



thank the OAN staff for assistance during the observations, Dr. J.P 
Sareyan for measuring the OAN dead-time corrections, and PAPIIT 
INI 10102 for the funds provided. Some support from the Fund for 
Scientific Research - Flanders (Belgium) through project G.0178.02 is 
also acknowledged. Part of these data were acquired with equipment 
purchased thanks to a research fund financed by the Belgian National 
Lottery (1999). Extensive use was made of the ADS and Simbad (op- 
erated by the Centre de Donnees Astronomiques, Strasbourg, France) 
data bases. 

References 

Antonello, E. & Mantegazza, L., 1982, IBSV 2152, 1 

Barnes, T. G, Fekel, E, & Moffet, T. J., 1977, PASP 89, 658 

Beech, M., 1989, Astroph. Space Sci. 152, 329 

Cowley, A. P., 1976, PASP 88, 95 

Cowley, A. P. & Bidelman, W. P., 1979, PASP 91, 83 

Cuypers, J. 1987, Medelingen van de Koninklijke Academie voor 

Wetenschappen, Letteren en Schone Kunsten van Belgie, Jg. 49, 

Nr. 3, ed. Paleis der Academien, Brussel, p. 21 
Danziger, I. J. & Dickens, R. J., 1967, IBSV 150, 1 
Eggen, O. J., 1979, ApJS 41, 413 
Elliot, J. E„ 1974, AJ79, 1082 

ESA 1997, The Hipparchos and Tycho Catalogues, SP 1200 (Vol. 1) 

Fekel, F. C. & Bopp, B. W., 1977, PASP 89, 216 

Fremat, Y., Lampens, P., Hensberge, H. et al., 2004a, In: Variable Stars 

in the Local Group, IAU Coll. 193, ASP Conf. Series 310, eds. D. 

Kurtz, K. Pollard, p.395 
Fremat, Y, Lampens, P., Hensberge, H., 2004b, MNRAS 356, 545 
Garrido, R., 2000, In: Delta Scuti and Related Stars, ASP Conf. Ser. 

210, eds. M. Breger & M.H. Montgomery, 67 
Grigahcene, A., Suarez, J.C., Dupret, M.-A. and Garrido, R., 2005, In: 

Tidal Evolution and Oscillations in Binary Stars, ASP Conf. Ser., 

eds. Claret, A., Gfmenez, A. & J.-P. Zahn (in press) 
Gr0nbech, B., Olsen, E.H., Stromgren, B., 1976, AAS 26, 155 
Hartkopf W.I., Mason B.D., McAlister HA. et al., 2000, Astron. J. 

119, 3084 

Hauck, B., 1973, in: Problems of Calibration of Absolute Magnitudes 

and Temperature of Stars, IAU Symp. No. 54, eds. Hauck, B. & 

Westerlund, B.E., p. 117 
Hauck, B., 1993, in: The MK Process at 50 Years, ASP Conf. Ser. 60, 

eds. Corbally, C. J., Gray, R. O., Garrison, R. E, p. 157 
Hoffleit, D. & Jaschek, C, 1982, The Bright Star Catalogue, Yale 

Univ. Obs. 
Jerzykiewicz, M., 1993, AAS 97, 421 

Kunzli, M., North, P., Kurucz, R.L. et al., 1997, AAS 122, 51 
Kuschnig, R., Weiss, WW., Gruber, R. et al., 1997, AA 328, 544 
Lampens, P. & Boffin, H. 2001, In: Delta Scuti and Related Stars, ASP 

Conf. Ser. 210, eds. M. Breger & M.H. Montgomery, 309 
Lampens, P., Garrido, R., Parrao, L. et al., 2004, In: The environment 

and evolution of double and multiple stars, IAU Coll. 191, Rev. 

Mex. de Astronorma y Astrofisica 21, eds. C. Allen & C. Scarfe, 

73 

Lampens, P., Fremat, Y, Uytterhoeven, K. et al., 2005, In: Tidal 
Evolution and Oscillations in Binary Stars, ASP Conf. Ser., eds. 
Claret, A., Gfmenez, A. & J.-P. Zahn (in press) 
Lester, J. B., Gray, R. O. & Kurucz, R.L., 1986, ApJSS 61, 509 
Maitzen, H. M., Pressberger, R. & Paunzen, E., 1998, AAS 128, 573 
Montgomery, M. & O'Donoghue, D., 1999, Delta Scuti Newsletter 
13, p28 

Mermilliod J.C. & Mermilliod M„ 1998, The General 
Catalogue of Photometric Data, new version 
I http://obswww. unige. cli/gcpd/gcpd. html) 



Lampens et al.: Photometric investigation of the S Scuti star DG Leo 



9 



Meylan, G., Python. M. & Hauck, B. 1980, AA 90, 83 

North, P., Jaschek, C, Hauck, B. et al., 1997, in: 'Hipparcos- Venice 

'97',ESA SP-402, p. 242 
Schoenaers, C. & Cuypers, J., 2004, In: Variable Stars in the Local 

Group, IAU Coll. 193, ASP Conf. Ser. 310, eds. D. Kurtz, K. 

Pollard, p. 283 
Sperl, M., 1998, Manual for PERIOD98 

I http://www. astro, univie. ac. at~dsn/dsn/PERlOD98 ) 
Schwarzenberg-Czerny, A., 1991, MNRAS 253, 198 
Rosvick, J. M. & Scarfe, C. D., 1991, PASP 103, 628 
Wilson, R.E., 1994, IAPPP Communications, spring 1994 



10 Lampens et al.: Photometric investigation of the 6 Scuti star DG Leo 

Table 2. Mean differences between variable/comparison 2 and comparison 1 star per season and per filter 



Data type 


Campaign 


Nr 


Filter u 


Filter v 


Filter b 


Filter y 










Mean (St.dev.) 


Mean (St.dev.) 


Mean (St.dev.) 


Mean (St.dev.) 


V- 


CI 


OAN 2002 


429 


-2.2251 (0.022) 


-2.2572 (0.014) 


-1.6962 (0.012) 


-1.3136 (0.011) 


C2 


-CI 


OAN 2002 


420 


-1.1529 (0.015) 


-0.7820 (0.006) 


-0.3107 (0.005) 


-0.0353 (0.004) 


V- 


CI 


OSN 2002 


748 


-2.2205 (0.011) 


-2.2608 (0.009) 


-1.7454 (0.008) 


-1.3279 (0.007) 


C2 


-CI 


OSN 2002 


750 


-1.0932 (0.009) 


-0.7351 (0.004) 


-0.3120 (0.003) 


-0.0424 (0.004) 


V- 


CI 


OAN 2003 


269 


-2.2216 (0.016) 


-2.2835 (0.011) 


-1.7666 (0.010) 


-1.3264 (0.008) 


C2 


- CI 


OAN 2003 


269 


-1.1408 (0.011) 


-0.7843 (0.002) 


-0.3155 (0.002) 


-0.0340 (0.002) 


V- 


CI 


OSN 2004 


144 


-2.2317 (0.010) 


-2.2697 (0.008) 


-1.7509 (0.007) 


-1.3336 (0.006) 


C2 


-CI 


OSN 2004 


144 


-1.1052 (0.011) 


-0.7407 (0.003) 


-0.3160 (0.004) 


-0.0459 (0.003) 



Table 7. Mean standard colours and global atmospheric parameters of DG Leo in two photometric systems 



a) Stromgren 


V 


b-y 


ml 


cl 


(b-y)o 


m 


Co 


P 


M v 


Teff 


logg 


[M/H] 




mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


mag 


K 


dex 


dex 




6.0788 


0.1462 


0.2210 


0.9328 


0.140 


0.223 


0.932 


2.788 


1.20 


7540 


3.59 


0.32 




(0.012) 


(0.003) 


(0.005) 


(0.008) 










(0.30) 


(100) 


(0.12) 


(0.12) 


b) Geneva 


V 


[U-B] 


[V-B] 


Bl-B 


B2-B 


Vl-B 


G-B 




M v 


Teff 


logg 


[M/H] 




mag 


mag 


mag 


mag 


mag 


mag 


mag 




mag 


K 


dex 


dex 



6.0817 1.6083 0.6467 0.9625 1.4083 1.3656 1.7780 
(0.012) (0.003) (0.006) (0.002) (0.004) (0.005) (0.006) 



1.08 7340 3.90 0.21 
(0.31) (100) (0.12) (0.12) 



